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Abstract 

Here is proposed the qualitative and quantitative solution of long-standing problem in astrophysics: 
the origin of 'knee' in Galactic cosmic rays (CRs). The paper deals with calculations based on the 
standard model of cosmic ray acceleration in Supernova remnants (SNR) and on the latest astronomical 
data concerning the variety in Supernovae types and energies as well as the sites of their explosions. 
The results obtained allow explanation of the main features of all particle cosmic ray spectrum: the 
'knee' around 3 PeV, the change of a slope by 0.3-0.4 in this point and the start of dip around 1000 
PeV. The results are based essentially on the premise that energy of accelerated CRs in the given SNR 
' is proportional to the kinetic energy of this explosion, that stresses the input of high energy explosions 

^ ■ and suppresses the contribution of low energy explosions to the total CR flux. 

m ; 

m ! 1 Introduction 

O ; 

Supernovae (SNe) represent the catastrophic explosions that mark the end of the life of some stars. It 
^■f is well known, that the mechanical energy input to the Galaxy from each supernova is about 10 51 erg. 

so with a rate of about 0.01-0.03 year -1 total power of SNs in our Galaxy is enough to provide the 
+3 . total energy of Galactic cosmic rays (GCR) ~ 10~ 12 erg/cm 3 (Berezinsky et al. p.990 ). It is shown 

that there exists a mechanism needed for the channeling about 10% (or even more (Berezhko et al. 
|2000aj ) of the mechanical energy of the explosion into relativistic particles. Considerable collective 
efforts have been made during the recent years to clarify the mechanism of CR acceleration in SRs 
(Drury et al. |2001j ). Theoretical progress is connected with the development of the kinetic nonlinear 
theory of diffusive shock acceleration (Berezhko et al. [T3Mj . pi??] . |2000aj . |2000b| : Berezhko& 
Ellison [T999] : Ellison et al. [TMT] . [20*00] : Drury et al. [207)1] : Malkov et al. fMT\\ ) and key advance 
has been made due to improved understanding of the nonlinear reaction effects on the shock structure. 
This theory is able to explain not only the main characteristics of the observed all particle cosmic ray 
spectrum up to an energy of 10 14 — 4 • 10 14 eV (Berezhko et al. {2QQQaj ; Drury et al. 2001]), but also 
heavy element abundances in CR flux relative to the solar system (Ellison et al. |1997j ). 
The standard model predicts (Drury et al. |2U01| 1: 

1) The power-like and approximately similar spectra of various nuclei of CRs beyond an energy of 
100 GeV/n. Spectra of heavy component can be slightly harder (A7 ~ 0.3) than proton one due to 
more effective acceleration of dust grains and ions (Ellison et al. [1997j ). 

2) The slope of power-like source spectrum is 7 so «r ~ 2.0 — 2.1. May be spectra have a 'curvature' 
being even harder before the maximal energy of acceleration, if nonlinear reaction effects are strong 
(Berezhko et al. j2()()()aj 1. (Ellison et al. [MI]). 
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3) The maximal energy of accelerated, particles E max (cut off energy) is ~ Z x 10 14 ev for the 
average SNe exploding into the average interstellar medium (ISM). 

4) There is a possibility to move the maximal energy to higher energies assuming nonordinary 
medium, for example, circumstellar medium (CSM) for any class of explosions: explosions into wind 
of Wolf-Rayet stars or explosions into superbubbles (Bykov & Toptygin 1997.), see also reviews 
(Ptuskin 2000 ; Biermann 2000 ). This effect is connected mainly with the higher magnetic field in 
the stellar or superbubble interior. 

5) The actual source spectra inferred from observations after propagation corrections is 7 SO w = 
jobs — A7. The value of A7=0.6 for simple diffusion model and A7=0.3 for model with reacceleration, 
A7=0.8 for model with Galactic wind (Jones et al. J2001J). So the value of A7 is not well known yet. 

The most 'nasty' problem, as it was called in (Drury et al. [2001 ). is the knee problem. The 
standard picture makes a clear prediction, that the GCR spectrum should start to cut off at energy 
about 10 14 eV or less for all species and drop exponentially as one goes to higher energies. Only 
some subclass of SNR can provide the knee particles while most SNRs have spectra cutting off at 
considerably lower energies (Reynolds & Keohane 1999 ). 

The upper limit of accelerations is determined essentially by the product of shock radius Rsk, 
shock velocity Vsk, ejected mass Mej, remnant age Tsnr, explosion energy E51, usually normalized 
to value of 10 51 erg. All these values are connected with each other and vary from explosion to 
explosion. E max can be expressed by the simple formula, if only Sedov phase of SNR expansion is 
considered (Ellison et al. |1997j ): 

M W U \AJ U-3/xGAciW Vl0 51 ergj [lO^km/sJ 1 h 

A week dependence on Vsk and on £"51 = E$n/W 51 erg, but strong dependence on magnetic field 
B of ISM can be seen there. This formula predicts E max ~ 100 • Z TeV for E51 = 1 and for the 
density of protons in medium uh ~ 0.3 cm -3 . An accurate determination of E max in frames of a more 
complete kinetic theory (Berezhko et al. }20Q0a| , |2000b ) taking into account timescale evolution of 
the shock in both Sedov phase and in free expansion phase of SNR evolution, gives, within a factor of 
2, the same results (Drury, |2001| ^. 

The main idea of this work (it seems very natural) is an attempt to get the average cosmic ray 
particle spectrum taking into account the distribution of SNRs in types, explosion energies, type of 
interstellar and circumstellar medium around it. In this case the total CR flux F(E) can be expressed 
by the formula 

Nz Ntp E 51 max 

F ( E ) = zZzZ j ^j(E 5 i)G(E,E max (E 51 ),B,Z)dE 51 , (2) 

where J2i 1S the summation by different cosmic ray nucleus groups Nz, J2j is t ne summation by 
different types of SNe explosions Ntp. G(E, E max (E51, B, Z)) is the average spectrum of comic rays 
in every explosion, which was approximated by the power low: 

G(E,E max ) = IoE-\ (3) 
7 = 2.0 in the interval 10 GeV < E < E max /5; 

7 = 1.70 in the interval E max /5 < E < E max ; 
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7 = 5 in the interval E > E % 



This spectrum shape takes into account nonlinear reaction of CRs on shock structure (Berezhko et 
al. |1997j . |2000a| ): the decreasing of 7 before E max . E max depends on E51, Vsk, B, Hh accordingly 
to formula (1). For every type of explosion the value of B and nn can be different. 

Intensity of CRs produced in every SNR is found from the condition 



As it was shown in (Berezhko |2U0Uap the portion of SNR kinetic energy transforrmed to CRs can 
reach even 30%. 

The most indefinite distribution is the distribution in .E51 for various types of SNe. The third 
section will be devoted to this problem. In the second section the latest astronomical observational 
data about variety in SNe explosions will be reviewed. In the fourth section we present the numerical 
results of calculations. In the fifth part of the paper we discuss physical interpretation of the results 
obtained. 

But we would like to note in advance that the diversity of SNe by explosion energies results in very 
unexpected and promising conclusion: the knee region moves to higher energy by a factor about of 
5-10, though we use only the standard model of CR acceleration and add to it the latest astronomical 
data on Supernovae explosions. 



Detailed observations of a growing number of Supernovae show that the nature of this phenomena 
is complex (Turatto 2003] ) . many new peculiar events discovered during last years display the wide 
range of luminosities, expansion velocities, chemical abundances, that is evidence for large variations 
in explosion energy and in the properties of their progenitors (Hamuy 2003 ). 

The thermonuclear explosions of accreting white dwarfs that explode as they approach the Chan- 
drasecar mass (~ 1.4M©) produce type la SNe, which owing to their high luminosity and accurate 
calibration are successfully used for the determining the geometry of Universe (Leibundgut 2000 ) , as 
a 'standard candles'. The energies of the explosions are practically fixed. 

Core collapse supernovae CCSNe (SNIb/c, SNII) are thought to be the gravitational collapse of 
massive stars (M > 8MQ), which makes the neutron stars as a compact remnants. CCSNe prove 
to comprise the most common general class of exploding stars in the Universe and they come in a 
great variety of flavors (Hamuy 2003 J. Even subclass of 'normal' SNII : plateau SNII-P and Linear 
SNII-L demonstrates a wide range in explosion energy, from 0.6 to 5.5 (xlO 51 ) erg among classical 
SNII (Hamuy, 2003). The ejected masses also are in broad range between 14 and 56 MQ (Hamuy 
2003] ). Despite the great diversity displaced by SNII-P, these objects show a tight luminosity- velocity 
correlation. This result suggests that, while the explosion energy increases so do the kinetic energies 
(Hamuy 200$]). These stars explode as isolated stars. 

A distinct class of SNIIdw can be identified which are believed to be strongly interacting with 
a 'dense wind' produced by SN progenitors prior explosion. These SNe have strong radio emission 
caused by the interaction with circumstellar medium (Chevalier |1D9S| ) (mass loss rate is ~ 10 _4 M©) 
for their progenitors. SNIIdw comprise ~ 30% of all CCSNe (Hamuy 2003 ). When the narrow line 
is present, the SN is classified as Iln ('narrow'). The strong degree of individuality is seen in their 
spectra, but despite the great photometric diversity among SNIIdw, these objects share the property 
of being generally more luminous than the classical SNII (Hamuy 2003 ). Among this type SNe one 




(4) 



2 Variety in Supernovae 
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event SN 1997cy is much more energetic than any other SNII (E ~ 30 • 10 51 erg (Turatto ^2000] ))• SN 
1997 cy and its twin SN 1999E (Rigon 2003 ) are associated to GRBs. As in the case of others these 
events show strong ejecta-CSM interaction with explosion energies as high as 3 ■ 10 52 ergs (Turatto 

mom). 

Hydrogen-deficient supernovae SNIb and SNIc are associated with the gravitational collapse of 
massive stars (may be Wolf-Raet stars), which have lost their hydrogen envelope during the phase of 
strong wind. In the case of SNIc most of the helium is gone as well. There is yet no direct observational 
proof for binary companions in SNIb/c, but this seems likely (Turatto 2002 ). 

In the past few year 3 SNe (1997ef, SN 1998bw, SN 2002 ap) have been found to display very 
particular spectra: they are extremely smooth and featureless, that can be interpreted as the result of 
unusual expansion velocities (Hamuy |2003j ). This fact suggests that these objects are hyper-energetic 
so they are called 'hypernovae'. The estimated energies of explosions are very high: 7 • 10 51 for 2002 
ap (Mazzaly [US]), 8 • 10 51 for 1997ef (Nomoto et al. fIM)\ ). 60 • 10 51 for 1998 bw (Nomoto et al. 
2000 ). The estimated expansion velocity of this object is as high as > 30000 km/s (Turatto et al. 
2002 ). SNe 1998bw was not only remarkable for its great expansion velocity and luminosity, but also 
because it exploded at nearly the same location and time as GRB980425 (Galama et al. 1998). 

M. Hamuy in the review (Hamuy j2003j ) makes the following conclusions. Despite the great 
diversity of core collapse SNe several regularities emerge which suggest that 1) there is a continuum 
in the properties of these objects, 2) the mass of the envelope is one of the driving parameters of the 
explosion, 3) the physics of the core and explosion mechanism of all core collapse SNe are not to be 
fundamentally different, regardless of the external appearance of the supernova. 

The great observational diversity of CCSNe has not been fully understood even if it clearly involves 
the progenitors masses and configurations at the time of explosion. Whereas SNII-P are thought to 
originate from isolated massive stars, a generalized scenario has been proposed in which common enve- 
lope evolution in massive binary systems with varying mass ratios and separations of the components 
can lead to various degree of stripping of the envelope (Nomoto et al. 1,1995 ). According to this 
scenario the sequence of types IIL lb Ic is ordered according to a decreasing mass of the envelope 
(Turatto et al. fKW\ ). 

3 How to estimate V(E 51 )? 

For this purpose we use the data (Richardson et al. [2)02 ), where a comparative study of the absolute- 
magnitude distributions of supernovae has been done. They used the Asiago Supernova Catalog (ASC), 
where by June, 2001 the number of events had increased to 1910, but number of events, for which 
this study is possible is much smaller. As far as the SNII types, either the number of events for 
which accurate peak apparent magnitude have been reported remains small (SNe lb, Ic, II-L, Iln) or 
the intrinsic dispersion in the peak absolute magnitude is large (SNII-P). For the absolute- magnitude 
distribution (in blue filter Mb), the authors consider only SNe within 1 Gpc. These distributions for 
different types of SNe are presented in Fig. 1. 

The analysis shows that (Richardson et al. [2002 ): 

1) At least 7 of 31 SNe in our Galaxy and in galaxies within 10 Mpc appear to have been sublu- 
minous (Mb> -15). Assuming that there is observational bias, it appears that more (perhaps much 
more) that 0.2 of all SNe are subluminous, but this fraction remains very uncertain. 

2) Only 20 of 297 extragalactic maximum-light SNe appear to be overluminous (Mb< -20), but 
it has become clear, that they exist. The absolute-magnitude dispersion of SNIbc has increased in 
comparison with the previous works due to the discovery of some rather luminous events. The SNe 
Iln are on average the most luminous type of core-collapse SNe. Considering the strong observational 
bias in favor of them, it is safe to conclude that the fraction of all SNe that are overluminous must be 
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Figure 1: Absolute-magnitude distributions of various types of SNe from (Richardson et al. 
[200|). 
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lower than 0.01. 

3) The authors have approximated absolute-magnitude distributions for every types of SNe by 
Gaussian functions with parameters, listed in Table 1. They consider also 'intrinsic' distributions 
obtained taking into account not only Galactic extinction, but also calculating extinction distributed 
for each SN type, averaged over all galaxy inclinations. Moreover, they divided Ibc into two luminosity 
groups: bright and normal. The II-L group was also divided into two groups. 

Table 1: Parameters of Gaussian distributions for 5 (7) main types of SNe from (Richardson |2002| ). 



SN Type 


< Mbim > 




Statistics 


Normal la 


-19.46 


0.56 


111(30) 


Total Ibc 


-18.04 


1.39 


18 


Bright Ibcl 


-20.26 


0.33 


5 


Normal Ibc2 


-17.61 


0.74 


13 


Total IIL 


-18.03 


0.90 


16 


Bright IIL1 


-19.27 


0.51 


4 


Normal IIL2 


-17.56 


0.38 


12 


Total IIP 


-17.00 


1.12 


29 


Total I In 


-19.15 


0.92 
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As was pointed out in (Hamuy 20031); ^ ne physics of the core and explosion mechanism of all 
core collapse SNe are not to be fundamentally different, so one can expect the presence of correlations 
between average absolute magnitude Mb for given SNe group and average energy of explosion for 
it. At least on the one hand the tail of most energetic explosions of Hypernovae is attributed to 
SNIbc, SNIIn (Hamuy 2003]) and on the other hand they are most luminous (see Table 1). Actually 
thermonuclear explosions of SNIa should be excluded. So we have made an attemt to estimate the 
average dependence Esi(Mb). 

For this purpose we use calculations of Nadyozhin (Nadyozhin 2003 ), where he makes predictions 
(on the base of hydrodynamical modelling of type II plateau supernovae light curves in the frames of 
the LN85 model (Litvinova & Nadyozhin |1983j )) for correlations between three observable parameters 
(the plateau duration At, the absolute magnitude My measured in V-filter and photospheric velocity 
Vph at the middle of the plateau) on the one side and three physical parameters (the explosion energy 
£?5i, the mass envelope expelled M e j and presupernova radius R) on the other side. According to 
LN85, the following approximate relations can be used to derive E51, M e j, R from observations. 

\ g E 51 = 0.135M y + 2.34 lg At + 3.13 lg Vph - 4.205 

lg Mej = 0.234M y + 2.91 lg At + 1.96 lg Vph - 1.829 (5) 

lgR = -0.572M y - 1.07]gAt-2.74]gVph- 3.350, 

where Mej and R are in solar units and Vph in 1000 km s _1 . The analysis of real 14 events SNe 
II-P shows (Nadyozhin 2003]): for this events the expelled mass, explosion energy and presupernova 
radius remain approximately within the limits Mej ~ 10 4 30MO, R ~ 2OO^-6OOi?0, E$\ ~ 0.6 4- 2.7. 
We rewrite these formulae to exclude the parameter At and obtain the following simple relations 

lg Esi = -0.43M y - 0.77 lg R + 0.52 lg Mej - 5.83 (6) 
lg Vph = +0.57 lg E 51 - 0.06 lg R - 0.48 lg Mej + 1.32 (7) 
For appropriate values Mej=20, R=200, we have 



l g E 51 = -0.43M y - 7 
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And we will use (7) further as a zero approximation. This is a key dependence for our calculations 
and the sensitivity of the results to this dependence will be discussed in fourth section. 

\gVph = +0.57 lg£ 5 i + 0.57 (8) 

In formula (1) the maximal energy of accelerated CR depends slowly on parameters of the explosion 
E§i and Vph as E max ~ 100(£ , 5i • Vph) 1 / 3 TeV. Using (7), (8), expressions we obtain E max (E^,i) 

lgE max = 0.52 • lgE 51 - 0.02 IgR - 0.16 lg Mej + 2.43 (9) 
For the same appropriate parameters Mej=20, i?=200 we have 

lg E max = 0.52 * lg E B1 + 2.1 (TeV) (9) 

So we obtain the dependence 

E max = 120£ 5 °i 52 TeV, (10) 

needed for our calculation by formula (2). 

For the estimation of ^f(E^i) we use dependence (7) and the sum of Gaussian distributions in 
Mb from Table 1 (Richardson al. 2002 ) for various types of SNe with average parameters listed 
in Table 1. Since lg-Esi depends liniary on Mb (7), the distribution in Mb can easily be converted 
to distribution in dN / d\gE^\. It will be also the sum of Gaussian functions with parameters: < 
\gE 51 >j= -0.43 < Mb >j -7, aj(lgE 51 ) = 0A3aj(Mb), where j means the type of SNe. Then 
dN/dlgE^i can be transformed into dN/dE^i = ^>(E§i) distribution. 

Here it is important to note that for SN la the physics of the explosion is absolutely different. So 
for SNIa we use < E§\ >~ 0.8 • 10 51 erg with dispersion of 20%. 

The final distribution ^f(E^i) for all types of SNe is presented in Fig. 2. As can be seen the most 
of events have energies (0.5 + 3-)10 51 erg. A remarkable peculiarity of the distribution is a very long 
tail toward the high energies up to I?5i=50. The reality of this is a key problem of our calculations 
and will be discussed in detail in the last section. 

Further we consider the most abundant sites for SNe. SN explosions are not random in the Galaxy, 
and all of them show strong spatial concentration to the centers of galaxies and to the arms in spiral 
galaxies (Bergh 1997 ). The degree of the concentration is larger for short-lived massive stars, initiated 
SNII, SNIbc explosions. SNe of type la occur in galaxies of all types, even far from the place of active 
star formation zones, due to the long evolutionary history in binary systems. Type Ibc appears only in 
spiral type galaxies and have been associated with a parent population of massive stars, perhaps more 
massive than SNell (Turatto 2003 ), and demonstrates even stronger concentration to the centers of 
galaxies, than other types (Bergh [1997] ^ . 

Numerous regions of very hot and rarefied gas with the temperature T = 5 ■ 10 6 K and proton 
density uh = 0.003 cm -3 occupy 50% of the volume of spiral arms of Galaxy (Kononovich & Moroz 
2001 ). A source of heating is thought to be the activity of young stars, first of all supernovae 
explosions. So for young massive stars this site of SNe can be selected as most probable for SNIbc, 
SNII. 

This variant of interstellar medium with parameters T = 10 6 K, nn = 0.003 cm" 3 , magnetic field 
B = 3yuG was used in (Berezhko |2QQQbJ ) being named 'hot phase' (Berezinsky et al. ^90 ). In this 
variant the maximal predicted energy is ~ Z ■ 300 TeV. 

Numerous regions of neutral gas HI can be divided into two parts (Kononovich & Moroz |2001j ^: 
the clouds of gas and dust with tlh = 10 cm -3 and T = 80K, that occupy relatively small volume 
1%, and intercloud regions, that occupy 50% of the volume of spiral arms with nu = 0.1 cm" 3 and 
T = 10 4 K. These latter regions were used in calculations (Berezhko j2000bj ) as a most probable site 
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Figure 2: The distribution ty(E 5 i) for sum of all types of SNe, converted from Mb-distribution, 
presented in Fig.l. 

of SNe with T = 10 4 K, nn = 0.3 cm -3 , B = 5/iG ('warm phase'). Variant with 'warm phase' gives 
the maximal energy of acceleration about Z ■ 100 TeV. We choose this site as most probable for SNIa. 

Besides that there exists the temporal correlations resulting from the concentration of the majority 
of core collapse SN progenitors into OB associations. The explosion of the first SN among such an 
association is to be followed by several tens of others. This results in the formation of a multiple 
supernova remnants, powered by both the SN explosions and the strong winds of Wolf-Rayet stars in 
the OB associations, which grows as a large bubble of hot, tenuous plasma known as a superbubble 
SBs(Tomisaka [1335] : Korpi [T333] b 

The SB acceleration model has been developed by (Bykov & Fleishman 1992 ), (Bykov & Uvarov 
1999]). They (Bykov & Toptygin 1997 ) estimated the maximal energy of accelerated nuclei as 10 18 
eV due to reacceleration effects, in the presence of magnetic field in the bubble interior of an order of 
30 /xG. In this model the spectrum beyond the knee is dominated by heavy nuclei. 

Since SNIIn explode in circumstellar medium (in accordance with a definition) and can be found 
in superbubbles, we choose much stronger magnetic field B = 30 — 45/uG, that increases the maximal 
accelerated energy by a factor of 10-15 for this type of explosions. 

In the next section we present numerical results of calculations of all particle cosmic ray spectrum, 
using all above mentioned dependences and parameters needed for formula (2). 

4 All particles cosmic ray spectrum, numerical results. 

We divided all nuclei of cosmic rays in 5 rough groups of p, He, (C, N, O), (Mg, Si, Ne), Fe with 
relative intensities 0.36, 0.25, 0.15, 0.13, 0.15 correspondingly. This chemical composition takes into 
account the fact that heavy components have slightly harder spectra beyond 1 TeV, than light ones 
(Shibata 1995 ), the contribution of heavy nuclei is increasing toward higher energies in comparison 
with 'normal composition', obtained around 1 TeV (Erlykin & Wolfendale |2001| ). The spectrum 
shapes are the same for all components (see formula (3)). 

Here it should be noted that we don't take into account propagation effects and present source 
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Figure 3: Total proton spectrum (in relative units) generated by 7 different types of SNe with 
parameters from Table. 1. The contribution of each SNe type is presented separately. 



spectra that might be easily converted to observable spectra in accordance with a standard model 
Isour = lobs — A7. The value of A7 equals to 0.54 (or 0.3 or 0.8) depending on propagation model 
(Jones et al. [200l] b 

Fig. 3 presents the total proton spectrum generated by 7 different types of SNe with parameters 
from Table 1, calculated by formulae (1),(2) with Mb — E^i conversion function (7), E max — E^\ depen- 
dence (9). Only a portion of SNIa was slightly decreased to 30% (Wheeler (ZHDSj ). The contribution 
of each SNe type is presented separately in Fig. 3. 

It can be seen from Fig. 3 that: 

1) The contributions of most energetic explosions are stressed significantly due to expression (4) 
in our calculations: the total number and total energy of accelerated cosmic rays is proportional to 
the total kinetic energy of the explosion. For example, only 5% of cosmic rays are originated by SNIa, 
while they comprise about 30 % of all SNe. 

2) The total intensity of CR is practically formed by contributions of SNelbcl (bright branch, see 
Table 1) and SNIIn. 

3) The location of the 'knee' is determined by maximal energy of accelerated CRs in the most 
energetic explosions. 

All obtained regularities can be understood distinctly, if analytical expression for average value of 
Emax is written using formula: E max ~ 300 x E®f TeV (the value 300 Tev corresponds to 'hot phase' 
of ISM). The statistical weight of E max should be proportional to the total number of accelerated CRs 
(~ £51 (4). 



-E51 max 

J *(E 51 ) • 300 • E° 5i 52 ■ E 51 dE 51 

< >= E5imm E , ima x = TeV. (11) 

/ tt(£ 5 i) • E 51 dE 51 ^ 51 " 

Emmin 

For a symmetric distribution function ^SlAE^x) with a small dispersion, as for the case of SNIa, the 
factor < E^ > / < E^ > is close to 1. But for a very asymmetric function, as in Fig. 2, this factor 
can be by many times larger. For example, if we choose power like function ^(Esi) = 0.44 • E^' 7 in 
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Figure 4: Spectra of different cosmic ray nuclei. 



interval E 51 = 1 100 and *(£ , 5 i) = 0.44 in interval E 51 = 0.1 4- 1, the value of < E^j 5 > / < E^ > 
is 5. (This power like shape of ty(E$i) will be discussed in section 5). 

As it can be seen from Fig. 3, the first 'knee' in the proton spectrum is located around 3 PeV, 
while for the most probable energy of explosion E pro f, = 1 maximal energy E max ~ 300 TeV (the 
variant of 'hot phase' of interstellar medium in (Berezhko 2000b|). 

The second 'knee' is formed by the contribution of SNIIn explosions, because they are also very 
energetic (see Table 1) and they explode into specific circumstellar medium (in our calculations B = 
45/iG for SNIIn). 

Fig. 4 presents the spectra of different cosmic ray nuclei, calculated in our model (Mg-Si-Ne group 
was omitted from the Fig. 4). We would like to remind that E max (Z) = Z ■ E max (p) (1). Every 
nuclear component also has two 'knees'. The maximal energy of accelerated CRs is determined by Fe 
nuclei, originated in SNIIn explosions and located around 10 18 eV. 

In Fig. 5 the contributions of various types of SN explosions to all particle spectrum are displayed. 

The change of a slope of power-like spectrum £7 beyond the 'knee' (in the interval 3 PeV to 26- 3 
PeV) is determined by a portion of Fe nuclei w(Fe) in chemical composition of CRs before the 'knee' 
in the case when the contribution of SNIIn is small enough: 

Ig(lMFe)) 

" lg26 y J 

For w(Fe) ~ 0.15 — 0.20 the change of the slope is close to £7 ~ 0.5. 

If the contribution of SNIIn is noticeable, then £7 should be less. Besides, as can be seen from 
Fig. 5., the more diversity in explosion energies (as in the case of SNIIP), the more smooth behaviour 
of the spectrum beyond the 'knee' is observed. 

The calculations reproduce on the whole the main features of all particle spectrum, measured in 
EAS experiments (Sommers ISQHj): the 'knee' around 3 PeV, the change of slope by £7 ~ 0.3 — 0.5 
beyond the 'knee', start of dip around 10 18 eV. We do not touch the region beyond 'ankle'. 

In Fig. 6 the average mass number <lnA> in the main variant of calculations (when the chemical 
composition of CRs for SNIIn explosions is the same as for others) is presented in comparison with 
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Figure 5: The contributions of various types of SN explosions into all particle spectrum. 



the data obtained in the KASCADE experiment (see rev. of (Sommers |20Ul| )b The main variant 
predicts less heavy chemical composition, than in the KASCADE experiment in the range 10 16 - 10 17 
eV. But in accordance with (Bykov & Toptygin [1997] ^) the supernovae explosions in superbubbles can 
generate the CRs enriched by heavy nuclei. So a variant when proton and helium components are 
absent in CRs generated in superbubbles is also presented in Fig. 6. The experimental dependence 
< InA > {E) is between the two variants. 



5 The physical interpretation of results 

The key problem in our calculations is a reality of the long tail in explosion energy distribution ^f(E^i) 
and sensitivity of the 'knee' location to this tail. To analyze the second item, we show in Fig. 8 the 
proton spectrum of CRs with different upper limits of integrating E^max in formula (2). In the main 
variant we used E^\max = 80. Fig. 7 shows that the point of knee location is determined by the SN 
explosions with £"51 ~ 30-7-60. At E§\max = 10 the 'knee' is around 700 TeV, that is higher than 300 
TeV due to equation (11): < E max >= 300 < E\^ >, but it is not enough to reproduce experimental 
data. 

One of the most interesting recent developments in the study of SNe is the discovery of some very 
energetic supernovae, whose kinetic energy exceeds 10 52 erg - Hypernova (Nomoto et al. 2002 ). They 
can be directly identifyed with the explosions determining cosmic ray 'knee' region not only by energy 
of explosions, but also by the type of core collapse SNe. In our calculations, just SNIbcl (bright) 
and SNIIn give the principal contribution to the formation of knee region (see Fig. 3, 5). Among 7 
possible Hypernovae 5 have been recognized as type Ic (1998bw, 1997ef,1997dq, 1999as, 2002p) and 
2 as type Iln (1997cy, 1999E) (Nomoto et al. 2002). The Hypernova branch might be interpreted as 
follows. Stars with M > 20 — 25M© form a black hole as a compact remnant; whether they become 
hypernovae or faint SNe may depend on the angular momentum in the collapsing core, which, in turn, 
depends on the stellar winds, metallicity, magnetic field and binarity. Hypernovae might have rapidly 
rotating cores owing possibly to the spiraling-in of a companion star in a binary system (Nomoto et 
al. [H!]). 
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Figure 6: The average mass number <lnA> in calculations in comparison with the data 
obtained in KASCADE experiment. 
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Figure 7: The dependence of 'knee' location on maximal energy of SN explosions. 
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Figure 8: The integral distributions W(> E51) for all mentioned in the text variants of calcu- 
lations. 

Since the physical picture of the Hypernovae explosion can differ from other core collapse SNe, 
we considered the variant, when in key dependence £51 (M6) (6) the parameters of ejected mass and 
radius of progenitor have been chosen different for SNIbcl, SNIIn (see Table 1) and other SNe. In 
the 'main variant' Mej=20, R=200. The variant with Mej=4, R=80 for SNIbcl and SNIIn, but 
Mej=10, i?=600 for other types from Table 1 is called 'Hypernovae variant' in Fig. 8, Fig. 9. 

Besides that, we consider the variant which was historically the first to be used in our calculations 
(the 'first variant'), and was based on a more sophisticated estimation of the dependence (6): lg£si = 
— 0.63Mb— 10.6. Here only five groups of SNe from Table 1 without division of SNIbc and SNIIL into 
bright and normal classes were considered. 

In ideal case the experimental distribution ^f(E^i) should be used. As a zero approximation to 
this case we consider of the tables from (Hamuy, 2003 ) (where the physical parameters for types II, 
Ibc, Ildw real supernovae are presented) and constructed the integral distribution W(> £51), since 
our calculation is sensitive mainly to the long tail in energies. This distribution (see Fig. 8) is called 
'experimental variant', while in reality it depends on the basic theoretical premises and moreover it 
can be slightly distorted by the selection bias. 

As it can be seen from Fig. 8, the distribution has a power-like shape with the slope of -0.79. It 
means that differential distribution is ~ £ -1 ' 79 . The portion of events with £51 > 1 was normalized 
to 0.6. In the region of £51 = 0.1 — 1 we choose N(E^i) = const. 

The integral distributions W(> £51) for the above mentioned variants of calculations are presented 
in Fig. 8. In Fig. 9 we show all particle CR spectrum calculated for every variant of W(> £51). 

The contribution of events with energy £51 ~ 50 comes to ~ 2% in all cases, but the contribution 
of events around £51 ~ 10 is by a factor of 2-3 larger than in 'experimental' distribution. But this 
difference occured to be not crucial for all particle CR spectra presented in Fig. 9, because the knee 
is mostly sensitive just to the energy of explosion £51 ~ 30 — 50, as shown in Fig. 7. 

So one can conclude, that the portion of events responsible for the formation of the knee in the CR 
spectrum comes to ~ 2 ± 1% meaning that with a usual SN-rate of about 0.01 year -1 the Hypernova 
rate is about 0.0002 year -1 . If the total power of 100 SNe with average energy £51 = 1 in our Galaxy is 
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E, TeV 

Figure 9: All particle spectra, calculated with different variants of W{> E 5i ) presented in 
Fig.8. 

enough to provide the total energy of Galactic cosmic rays ~ 10 -12 erg/cm 3 (Berezinsky et al ^)90 ), 
then 2 SNe with E51 = 50 or 3 SNe with E§\ = 30 can also provide the total energy of Galactic cosmic 
rays. 

If the time of CR life in our Galaxy is about 3 • 10 year and Hypervae rate is 0.0002 year -1 , it 
means that about 6 • 10 3 explosions provide the CR intensity in our Galaxy. But for high energy CRs 
(around 1 PeV) the time of CR life T esc is much smaller due to decreasing of escape lenght X esc as 
E~ a (a = 0.54) beyond the energy 5 GeV (Jones et al. 2001 ). The number of explosions giving the 
dominating input to CRs around the E ~ 1 PeV might be rather small (~ 10 — 15) in hole Galaxy 
and only several explosions in the nearby region. 

During the last millennium 5 SNe are known to have exploded (1006, 1054, 1572, 1054, 1679(Cas 
A)), see review of (Weeler [2DQ3") ) . The supernova of 1006 was the brightest ever recorded, brighter 
than Venus and perhaps as bright as a quarter Moon. Of these events, only SN 1054 produced an 
observable compact remnant, the Crab pulsar. Although there is still some controversy, the events of 
1006, 1572, 1604 are generally thought not to have left any compact remnant, but to have resulted in 
complete disruption of the star (Weeler ,2003 ). 

6 Conclusions 

1. The main idea of this work to take into account the distribution of Supernovae in types and 
explosion energies results in a very promising conclusion: the maximal energy of accelerated 
CRs moves to higher energies by a factor of 5-10 in comparison with the value 100-300 TeV, 
predicted by standard model of CR acceleration, due to diversity in SN explosion energies. This 
conclusion is based essentially on the premise that the energy of accelerated CRs is proportional 
to the energy of an explosion, that stresses the input of high energy explosions and suppresses 
the contribution of low energy explosions to the total CR flux. It is worth to note here that 
rather significant diversity in energies (around one order) exists even inside the SNIIP group 
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(Nadyozhin [2003 ; Hamuy fKM\ ). 

2. The knee in CR spectrum can quantitatively be explained by the dominating contribution of 
SNIbcl, SNIIn explosions with energy of (~ 30 — 50) • 10 51 erg, that might be identified with 
Hyper novae. 

3. The estimated rate of these energetic explosions is about 1-3 per 10 4 years being enough to 
provide the total power of Galactic cosmic rays. As a whole about 6000 stars provide the CR 
flux at lower energies, but may be only several explosions give the dominating contribution to 
the CR flux around the knee region (nearby of local Solar system) . 

The latter conclusion may intersect with the idea, proposed by Erlykin and Wolfendale (Erlykin &: 
Wolfendale [1997] ). that a single nearby local SNR accelerates the particles and gives the dominating 
input (mainly by O and Fe nuclei) to the knee region. But the CR propagation in the Galaxy is 
absent for this SNR. In our model the most energetic explosions give the dominating contribution to 
the whole energy range of CRs and their propagation in Galaxy should be taken into account, while 
we do not consider here these effects. Fractal diffusion seems likely to be very usefull to help CRs to 
reach the Solar system (Lagutin et al. |2001| ). 
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